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Abstract 

We consider the chargino contribution to the CP violation in B — > <pKs decay in the minimal 
supersymmetric standard model at large tan/3. It is shown that the Wilson coefficient C$g of 
the chromomagnetic penguin operator can be significantly enhanced by the chargino-mediated 
diagrams while satisfying other direct /indirect experimental constraints. The enhanced Cs g allows 
large deviation in the CP asymmetry from the standard model prediction, especially it can explain 
the apparent anomaly reported by BaBar and Belle. 
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The measurement of CP asymmetries at B-factories is a powerful probe of new physics 
(NP). In the standard model (SM), the origin of CP violation is the phase <5ckm in the 
3x3 Cabibbo-Kobayashi-Maskawa (CKM) matrix of quark mixing. In general, there can 
be many new sources of CP violation in new physics models beyond the SM. 

The time-dependent CP asymmetry in the neutral B decays to CP eigenstates B — * /cp 
gives information on the two classes of CP violation C/ and Sf [1]: 

r(B(t) - fcp) - r(B(t) - / C p) 



A C p(t) 



r(B(t) -> /cp) + r(B(t) -> /cp) 

-C/ cos(Am B t) + 5/ sin(Am B t), (1) 



where Am# is the mass difference of the B system and B(B)(t) is the state at time t which 
was B(B) at t — 0. The CP asymmetries C/ and Sf are determined by 

Acp^-W (2) 

where (3(6d) is the contribution of SM (NP) to the phase in the B — B mixing and A(A) is 
the decay amplitude for B(B) — > /cp. 

The time-dependent CP asymmetry in _B — > J/ipK decay measured at BaBar [2] and 
Belle [3] shows the existence of CP violation in the S-meson system, and the current world 
average [1] 

sin 2(3 j/iPKs = S J/i)Ks = 0.734 ± 0.054. (3) 

is fully consistent with the CKM picture of CP violation of the SM. In the SM, the CP 
asymmetry in B — > 4>K S also measures the same angle angle (5 because the CP violation in 
the decay amplitude is suppressed by the Cabibbo angle A (f» 0.22): S^Ks = Sj/^k s + 0(\ 2 ). 

Recently, BaBar [4] and Belle [5] have announced the first measurement of time-dependent 
CP asymmetries in B — > <pKs- The weighted average of the mixing-induced CP asymmetry 

sin 2/V s = S+Kg = -0.39 ± 0.41 (4) 

differs from the SM prediction of (3) by 2.7a. The Belle collaboration has also reported the 
direct CP asymmetry 

C^k = 0.56 ± 0.43, (5) 
which is consistent with zero as predicted in the SM. 



Contrary to the B — > J/^K S decay where A/ A = 1 to a good approximation even in the 
presence of NP, the loop-induced process B — > 4>K S generally allows large deviation from 1 
in phase and modulus of A/A. There are already many works which explain the apparent 
deviation (4) in various NP models [6-10]. 

The measured sin2/3's from B — > rj'Kg, B — > (K + K~)K$ [5] and sin 2a from B — > 
7r + 7r~ [11] are consistent with the SM although the errors are large. We note that, however, 
B — > 4>K S is unique in that it doesn't have the tree-level amplitude in the SM unlike 
most other 5-decays. Therefore, it is not unlikely that the NP manifests itself only in the 
B -> <pK s decay [12]. 

The minimal supersymmetric(SUSY) standard model (MSSM) has many new CP violat- 
ing phases besides the CKM phase of the SM. With CP violating phases of order one the 
electric dipole moments (EDMs) easily exceed the experimental upper bounds by several 
orders of magnitude. In addition, the general structure of the sfermion mass matrices in 
the generation space leads unacceptable flavor changing neutral currents (FCNC) by gluino 
mediation. These SUSY CP and SUSY FCNC problems strongly constrain the MSSM pa- 



in the MSSM, it has been shown that (4) can be accomodated if there exist new flavor 
structures in the up- or down-type squark mass matrices [6, 10, 13]. In this paper, we will 
show that the chargino contribution can generate large deviation in S^k at large tan f3, even 
if CKM is the only source of flavor mixing. 

Specifically, we adopt a decoupling scenario where the masses of the first two generation 
scalar fermions are very heavy (> O(10 TeV)), so that the SUSY FCNC and SUSY CP 
problems are solved without a naturalness problem [14]. We also assume the flavor-changing 
off-diagonal elements of scalar fermions are vanishing to guarantee the absence of the gluino 
mediated FCNC. In this case the CKM matrix is the only source of flavor mixing, while 
there are new CP violating parameters /i, M 2 in the chargino matrix 



rameters. 





and A t in the scalar top mass-squared matrix 




m 2 ~+rrit+D L m t {A\ — //cot/3) 
m t (A t - n* cot (3) m| + m\ + D R 
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where D L = (1/2 - 2/3 sin 2 6 W ) cos2/3m| and D R = 2/3 sin 2 9 W cos 2/3m|. 

In this decoupling scenario, it has been shown that the light stop and chargino contribu- 
tions still allows large direct CP asymmetry in the radiative B decay up to ±16% [15]. It 
should be noted that the new contribution to B — B mixing is very small and this model is 
naturally consistent with (3) [15]. 

The decay B — > <pKs is described by the AB = 1 effective Hamiltonian [16]. The chargino 
contribution to the Wilson coefficients of QCD penguin operator is given by 

Cf = E ^(xf)M) k sP2(x Ik \ (8) 

I,k=l xf 

where x t u = m t/ m2 H±, x ik = mt±/m2 and P2(x) is the loop- function [12]. (xi^)kq is the 
stop(/c)-chargino(/)-down-quark(g Z/ ( i? )) coupling: 

(x/% = -vas^ + v;^— -~ — 



y/2mw sin^ ' 



(xf)k q - UnSf t —= - -, 

\/2m w cos (3 



(9) 



where U, V (S) diagonalize(s) the chargino (stop) mass matrix. The other Wilson cofficients 
of QCD penguin operators are simply related to C3 by C5 = C3, C4 = Cq = SC3, where we 
have neglected the box diagram contributions which are suppressed by A. In the analysis 
we decouple the charged Higgs contribution by assuming m H ± = 1 TeV. The effect of H ± 
will be mentioned below. The expressions for the Wilson coefficients of magnetic operators 
C* l8g can be found, for exmample, in [12, 17]. 

We note that the has the enhancement factors m~± / by the chirality flip inside 

the loop and can dominate the SM contribution. On the other hand the C3 ... 6 preserve 
the chirality and don't have such enhancement factors. In addition, due to the super-GIM 
mechanism, the chargino contribution to C 3i ... i6 are much smaller than the SM values. The 
contribution to chirality flipped operators C^ 8g ^ are suppressed by m s /m b . We neglect the 
contribution of electroweak penguin operators which are also expected to be negligible in 
our scenario. Therefore the large deviation in A/ A should be generated soley by C$ g in this 
scenario. 

To calculte the hadronic matrix elements we use the QCD factorization method in ref . [16] . 
In this approach it has been demonstrated that the strong phases of QCD penguin operators 
cancel out in the SM [7]. Since the NP contribution to C 3; ... j6 are negligible in our scenario, 



-3-2-10 1 2 3 

arg(C 8g ) 



FIG. 1: S<j>K as a function of avg(Cs g ) for \Cs g \ =0.33 (short dashed line), 0.65 (long dashed line) 
and 1.0 (solid line). 

the strong phase is small as in the SM. The required C§ (m b ) (C^(m b ) —0.147) for large 
deviation of CP asymmetries can be estimated from the approximate numerical expression 
for A [7, 16]: 



In Fig. 1, we show S^k as a function of arg(C 8s ) for |C 8g | = 0.33, 0.65 and 1.0. From this 
figure, we can see |C 8g | ~ 0.33 — 0.65 with large positive phase can accomodate the deviation 
within la. 

Since the chargino contributions to CV 7 and C& g have similar structures, one can think 
that large deviation in C$ g may result in the large deviation CV 7 . Too large deviation in CV 7 
will violate the measurement of B(B — > X s ^), for which we take 2 x 10~ 4 < B(B — > X s ^) < 
4.5 x 10~ 4 [18], because it is already consistent with the SM prediction. 

Due to the different loop functions, however, the two Wilson coefficients are not strongly 
correlated, and it is possible to have large deviation in C Sg while keeping \C 7j \ constrained 
to satisfy B{B — > -XVy). Since x R i n (9) is proportional to tan/5 for large tan/3, we need 
relatively large tan (3 to have sizable effects. 

The neutral Higgs boson h° which is lighter than the Z-boson at tree level has large 
radiative corrections [19] which can be approximated at large tan/3 as 



A K Yl V P*s V pb( a z + < + a s) 



p=u,c 



w -3.9 x 10" 4 (3.7e a21i + A.5C 8g ). 



(10) 




(11) 
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FIG. 2: Distribution of Cg values for tan/3 = 35 (a) and tan/3 = 60 (b) in the complex plane. See 
the text for other parameters. The points with S^k < are marked with + symbol. 

where v ~ 246 GeV. LEP II sets the lower limit on the SM-like Higgs boson mass: m^o > 
114.3 (GeV) at 95% CL [20]. This limit gives very strong constraint on the stop masses. 

At large tan/3(~ O(50)), the SUSY QCD and SUSY electroweak correction to the non- 
holomorphic couplings H*D C Q become very important. This gives large correction on the 
down-type quark masses and some CKM matrix elements [21] 



m h = 



Vjj = V 



eff 
Jl 



1 + 63 tan /3 
1 + eo tan (3 



(12) 



1 + 63 tan/3' 

where ? 3 « e + e Y y 2 t and (J I) = (13)(23)(31)(32). m b , Vjf are 6-quark mass and CKM 
elements measured at experiments, respectively. At one-loop and SU {2) L x U(1)y symmetric 
limit, the e and ey are given by 

2a, 



^Re' " 

3n \ m~ f 



where 



m%/\mg\ 2 , 



m 2 ~/\m 



91 ' 



(13) 

m~/\n\ 2 . The loop-function 



is given by j(x, y) = (j(x) — j(y))/(x — y) with j(x) — x\ogx/(x — 1). Note that these SUSY 
threshold corrections are not easily decoupled even for very heavy superparticles. 

For the numerical analysis we fix the CKM matrix by using \V US \ = 0.2196, \V c b\ = 
4.12 x 10~ 2 and |14&| = 3.6 x 10 -3 [20], leaving only #ckm as free parameter. For example, 
for 5 C km = tt/3, we get Am Bd = 0.491 ps _1 and sin 2/3 = 0.729, which are close to the 
experimental central values [20]. The free SUSY parameters in our scenario are tan/3, M 2 , 
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FIG. 3 



Correlation between C^k and S^k f° r tan (3 = 60. 



(j,, wig, mi and A t {m H ± = 1 TeV) (also mg,m^ are relevant for large tan/3). In Fig. 2, we 
show the distribution of C 8g in the complex plane for <5ckm = 7r/3, tan [3 = 35(60), m H ± = 1 
TeV, itlq = 0.5 TeV, rrig = 1 TeV, and =0.5 TeV. We scanned the other parameters 
as follows: 



We have fixed the phase on the gluino mass parameter such that arg(m^) + arg(yu) = ix 
to maximize the SUSY QCD correction. When scanning, we imposed the B(B — > X s ^) 
constraint and the direct search bounds on the (s)particle masses [20]: m h o > 114.3 GeV 
and m~±,m^ > 100 GeV. 

From Fig. 2, we can see that our scenario can easily accommodate the discrepancy (4). 
As mentioned above, we have chosen large value for the charged Higgs mass = 1 TeV to 
safely suppress the Barr-Zee type two-loop EDM constraints which are significant if tan (3 [22] 
is large and the pseudo-scalar Higgs boson (A ) is relatively light 1 . We have checked that 
for the smaller m#± the larger deviation in S^k is possible due to the cancellation in the 
real part with the chargino contribution [12]. Therefore Fig. 2 is a rather conservative result 
for S^k- We have checked that the allowed range for S^k is not sensitive to the change in 

1 The SUSY contribution to (g — 2) M is also small in this case. [23]. 



< m~ t < 1 TeV, < |//| < 1 TeV, 



< \A t \ < 2 TeV, < |M 2 | < 1 TeV, 



7T < arg(/x),arg(A),arg(M 2 ) < ir. 



(14) 



7 



£ckm if we impose the constraint (3). 

C^k is correlated with S^k- In Fig 3, the correlation is shown for tan/5 = 60. We can 
see that it can also accommodate (5) although it favors small negative C^k, which will be 
clarified by future experiments. 

The direct CP asymmetry in B — > X s 7 decay is also expected to be large when S^k has 
large deviation from the SM expectations. We have checked that actually it can be very 
large but it is not correlated with S^k- It is because S^k is determined by the complex 
phase on C Sg while A C p{B — > X s ^) is mainly controlled by that on CV 7 [24]. Because we 
have phase on three independent parameters //, A t and M 2 , the phases of CV 7 and C% g need 
not have strong correlations. 

The B(B — > <f>K) varies moderately over the parameter space we considered and seldom 
exceeds 15 x 1CT 6 , which is acceptable compared with the experimental measurements [4, 5]. 
Also the mass difference Am s in the B s — B s system is close to the SM expectation Am s ~ 
14.5 ps -1 in the most region of parameter space we considered, which may distinguish our 
scenario from other scenarios in refs. [10, 13]. 

For large tan/5 in the MSSM, through the Higgs-mediated FCNC the B(B S — > 
can be enhanced by a few orders of magnitude over the SM prediction [25]. Observation 
of this leptonic decay mode, for example at Tevatron Run II, would be clear signal of NP. 
However, this is possible only for relatively light A . Since the large deviation in the S^k, 
although it needs new CP violating phase (s) of 0(1), does not necessarily require light A 
as we have shown, these two decay modes can be complementary to each other in searching 
for the MSSM at large tan (5. 

In our scenario the sole source of large deviation in S^k is Cg g . It can have simultaneous 
effects on other decays such as B — > </></>, B — > nn, B — > Kir, etc [12]. 

In conclusion we have considered the chargino contribution to the CP asymmetries S^k 
and C^k of B — > 4>K S decay in the CP violating MSSM scenario at large tan/5. We 
have shown that through the enhanced Wilson coefficient C$ g of chromo-magnetic penguin 
operator by the large SUSY threshold corrections [21], there can be large deviations in the 
CP asymmetries. 
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